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Abstract 

We present p53-MDM2-Glucose model to study spatio-temporal properties of the system induced by 
glucose. The variation in glucose concentration level triggers the system at different states, namely, 
oscillation death (stabilized), sustain and damped oscillations which correspond to various cellular states. 
The transition of these states induced by glucose is phase transition like behaviour. We also found 
that the intrinsic noise in stochastic system helps the system to stabilize more effectively. Further, the 
amplitude of p53 dynamics with the variation of glucose concentration level follows power law behaviour, 
A s (k) ~ A: 7 , where, 7 is a constant. 
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Introduction 

Oscillations are inherent and inbuilt in living systems due to various fundamental molecular processes 
and coordinate basic biological functions and their mechanisms [T| to self-organize the complicated life 
processes [2] . The origin of these oscillations could be due to various control or feedback mechanisms and 
set of non-linearities which describe complicated non-linear activities in the system [3;. p53 is one of most 
important proteins in cellular system which contribute to the maintenance of the genomic integrity [4j, 
involves in various cellular activities, such as, cell cycle arrest, DNA repair, apoptosis and other cellular 
functions EH3, and exhibits oscillatory behaviour by interacting MDM 2 via feedback mechanism [7j. 
Regulating p53, variety of cellular stresses, for example, global DNA damage that cause abnormal or 
cancerous cells and repairing of stress-induced DNA damage, can be controlled [5118] . Even though, p53 is 
functionally inhibited in normal cells |6|9], the MDM 2, which is a negative regulator of p53, can activate 
p53 to induce stress in the cell mmm- 

The oscillatory behaviour exhibited by p53 due to negative feedback mechanism with MDM 2 can be 
regulated and induce stress by other molecular activities, such as, Ca 2+ , NO [12], MTBP nsj etc. and 
excess stress may lead to apoptosis mmm- However, if the stress induced by various stress inducers 
(NO, reactive oxygen synthase (ROS), ARF etc) is weak, on removing stress p53 activation may come 
back to its normal functioning m which can be done by MDM 2 by enhancing the ubiquitination of 
p53, as a result of which p53 level degrades jTTUTT)] . 

However, there are many open issues regarding the role of glucose in cellular activities, for example, 
switching of normal to stress states via glucose level in the cell, the way how excess stress caused by 
glucose lead the cell to apoptosis, the possibilities to estimate the critical concentration of glucose that 
cause the state transitions, the way how glucose effects p53 dynamics etc. We, in this work, study 
the impact of glucose on ROS activation which cause DNA damage in p53-MDM2-Glucose network to 
understand different states in the system and their transitions. The work is organized as follows. We 
describe p53-MDM2-Glucose model and its molecular interaction with numerical techniques in section 2. 
The numerical results are presented with discussions in section 3 and some conclusions are drawn based 
on the results we obtained in section 4. 
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Mathematical model of p53 — MDM 2 — Glucose network 

The model we consider is the extension of p53 — MDM 2 network m induced by two new and impor¬ 
tant molecules namely stress inducer glucose and ROS, which damages DNA (Fig. 1). In the model, 
p53 acts as a transcription factor and helps in the transcription of MDM2.mRNA via MDM 2 gene. 
MDM2_mRNA then synthesizes MDM2 protein through the translation process. p53 interacts with 
MDM2 by enhancing its degradation through its E3 ubiquitin ligase activity m and maintains low con¬ 
centration level at normal condition [TO] . We then consider glucose metabolism inside cell which leads 
to production of ROS JT5] and high ROS concentration triggers the DNA damage mm inducing stress 
to the system [2TJ[22l. The DNA damage activates ARF and this activated ARF interacts with MDM2 
forming ARFJMDM2 complex J23][23]. This activity of ARF blocks the MDM2 E3 ubiquitin ligase 
activity which promotes MDM2 degradation [23]. Our model (Fig. 1) consists of nine molecular species 
listed in Table 1 which undergo the reaction channels listed in Table 2. 

Table 1 List of molecular species 


S.No 

Molecular Species 

Description 

Notation 

1. 

p53 

Unbound p53 protein 

Xi 

2. 

MDM2 

Unbound MDM2 protein 

X2 

3. 

p53 MDM2 

p53/MDM2 complex 

x 3 

4. 

MDM2 mRN A 

MDM2 messenger RNA 

X 4 

5. 

Glucose 

Unbound Glucose 

x 5 

6. 

ROS 

Unbound ROS 

X6 

7. 

Dam DN A 

Damage DNA 

X 7 

8. 

ARF 

ARF protein 

Xg 

9. 

ARF.M DM2 

ARF/MDM2 complex 

Xg 


Cellular processes are basically noise induced stochastic processes [25] and this noise could be intrinsic 
due to random molecular interaction in the system |26| as well as extrinsic due to fluctuations of physical 
variables surrounding the system m- The intrinsic molecular interactions in our model system undergo 
seventeen reaction channels (Table 2) and can be well explained in stochastic manner [261I28U29] as follows. 

Consider a configurational state of our model system at any instant of time t is defined by a state 
vector, X(t) = [Ai(t), X 2 (t),..., Ajv(t)] T , where, {W} is the set of variables corresponding to populations 
of the molecular species (Table 1), N = 9, and T is the transpose of the vector. The time evolution of 
the configurational probability P(X;t) to have a transition from one configurational state X to another 
state X(t)' during the time interval [t,t + dt\ is given by the following Master equation, 

dP{X d f U) = ^P(A(t),t)T^,-^P(A',t)r^ (1) 

{X} {A'} 

where, {T} is the set of transition rates from one state to another. Since solving equation (JT]) is very 
difficult for our system, we follow Gillespie [28] to simplify it to Chemical Langevin equations (CLE) as 
follows. The number of reactions fired during the time interval [t, t +At] is a variable B{a) which depends 
on propensity functions (a) of the reactions, and one can impose two important realistic approximations 
in the large population limit to arrive at CLE. First in the limit, At —> 0, the values of a will remain 
constant during [t,t + At], which allows B to be replaced by statistically independent Poisson random 
variable. The second approximation is At —> oo which allows to approximate Poisson random variable 
by a normal variable, G with the same mean and variance. Now linearising G, we reach the following 



















3 


Table 2 List of biochemical reaction, Kinetic Law and their rate constant 


S.No 

Reaction 

channel 

Description 

Kinetic 

Law 

Values of rate 

constant 

References 

1 

fci, . 

X 4 - > X 4 + X 2 

MDM2 translation 

k\x 3 

4.95x 10 4 sec 1 

mmm- 

2 

k 2 

Xi - > Xi + X 4 

Synthesis of 

MDM2 mRNA 

k 2 x i 

1.0 x 10 ” 4 sec _1 

mmm- 

3 

X 4 - > <p 

Degradation of 

MDM2 mRN A 

k 3 x 3 

1.0 x 10 “ 4 sec _1 

Bainni. 

4 

ki X 

X 2 - > <p 

Degradation of 

MDM2 

k4X2 

4.33 x 

10 4 sec —1 

mmm- 

5 

<p - > Xi 

Synthesis of p53 

k 3 

0.78 molsec -1 

HEEZi. 

6 

k 6 

x 3 —> X 2 

Decay of p53 

k 6 x 3 

8.25 x 10 4 sec 1 

mmm- 

7 

1 fe7 v 

Xl + X 2 - > X 3 

Synthesis of 

p53_MDM2 com¬ 
plex 

k7X\X2 

11.55 x 

10 ~ 4 mol~ 1 sec ~ 1 

muim. 

8 

kg 

X 3 - > Xi + X2 

Dissociation of 

p53_MDM2 com¬ 
plex 

k s x 3 

11.55 x 

10 _ 6 sec _1 

mmm- 

9 

X fc 9, 

<p - > X 5 

Creation of Glucose 

kg 

k [mol~ l sec^ 1 ] 

Assuming the con¬ 
centration level of 
glucose vary within 
the cell. 

10 

&10, 

Synthesis of ROS 

^10*^5 

2 x 

10 ~ 3 mol~ 1 sec ~ 1 

Assuming ros pro¬ 
duction due to the 
glucose molecule. 

11 

fell 

Xq - > X 7 

DNA damage 

knxe 

5 x 10 4 sec 1 

D2II32IE3- 

12 

X 7 —^ <j> 

DNA repair 

ki 2 x 7 

2 x 10 5 sec 1 

mm- 

13 

&13, 

x 7 —> x 8 

ARF activation 

k 13 x 7 

3.3 x 10~ 5 sec 1 

mm • 

14 

X 2 + X s Xg 

Synthesis of 

ARF/MDM2 

complex 

ki4X2X 3 

1 

10 ~ 2 mol~ 1 sec ~ 1 

mm ■ 

15 

ki5. 

Xg - > X 8 

ARF dependent 

MDM2 degradation 

ki 5 Xg 

1 x 10 _ 3 sec _1 

0H23. 

16 

&16 v / 
x 8 —> (p 

Degradation of 

ARF 

kiexs 

1 x 10 _ 3 sec _1 

PEI. 

17 

fcl7 v / 

X 5 - > 0 

Degradation of 

Glucose 

knx 5 

1 x 10 _ 4 sec _1 

Due to half life of 
the molecule. 
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CLE for our model, 

dx i 
dt 
dx 2 
dt 


dx 3 
dt 


dx 4 
dt 
dx 5 
dt 
dx 6 

da: 7 
di 
dx 8 
dt 


dxg 

dt 


k 5 - krX\X2 + fc 8 X3 + 




Vfcfc - \JkiX\X2i2 + Vk^xliz 


k\x 4 — fc 4 x 2 + k 3 x 3 — £ 7 X 1 X 2 + fcs^ — 

H—4= \/fcrX4^4 - \fkix 2 ib + Vk%x^& - s/hjxix^^ 

W . 

+ Vh X 3 £ 8 - \Zk 14 X 2 Xs^9 


—k 6 x 3 + £ 7 X 4 X 2 — k s x 3 
1 


'Vv[ 


- \/kex 3 ^io + \Zk 7 xix 2 £,ii - y/k 3 x 3 £, 12 


k 2 x i — k 3 x 4 H— -j= 

W 


vfea?i£i3 — \/ k 3 X 4 t;i4 


kg — £10^5 — £17X5 + — j= 


y/kgi is — y/kigx^ie - 1/^17X5^17 


fci 0 x 5 - £nX 6 + 
fcnX 6 — £ 12 X 7 + 


1 

Vv 

1 

Vv 


VkigX5^is — \/k\ rxe/ig 


\/£i 1 X 6£20 - \A 12 X 762 

1 


k\ 3 X 7 - fci 4 X 2 X 8 + fci 5 Xg - £ i 6 X 8 + 




\/ kl 3 X 7 ^22 


— V kl4,X2X8&3 + 1/^15X9^24 + \/ ^16^8^25 


fcl 4 X 2 X 8 - fci 5 Xg + - 7 = 

w 


\ZVaX2X8^26 — 1 /^ 15^967 


( 2 ) 


(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 


(9) 

( 10 ) 


where, {£j} are random parameters which is given by, = limdt^oNi(0, 1)/Vdt and satisfy £i(t)£j(t') = 
6ij6(t — t'). When V —> 00 or {£} —>■ 0 the set of CLEs recover deterministic equations. 

The dynamical behaviour of the system can be well studied by simulating the set of CLEs given 
by equations (2)-(10) or corresponding deterministic differential equations by using standard 4th or¬ 
der Runge-Kutta method El for numerical integration of the set of stochastic or ordinary differential 
equations. 

We use the stochastic simulation algorithm due to Gillespie (26[ to calculate the time evolution of 
state vector of the system by simulating the reaction sets given in Table 2 with the parameter values. The 
algorithm is based on proposed joint probability density function P(r, fi) = T{t)R{^) to allow a transition 
to occur with time increment r (= [0, 00 ]) by finding the probability densities, T(r) at which a particular 
reaction is fired and R{V) which identifies the particular reaction /r = [ 1 , 2 fired at that time. 
From this hypothesis, r and /j, can be estimated computationally from the relations r = — ln[T(T)\/a 0 

and = uJoR(V) by defining two random numbers rq and r 2 for T(r) and R(V) respectively and here 
\-~\M 

a ° — Z^i=1 a i- 













































Figure 1. The schematic diagram of p53-MDM2-Glucose interaction network which involves 
feedback loops induced by ARF protein and regulated indirectly by Glucose. 
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Figure 2. Plots of a comparative behaviour of p53 and Ad DM2 as a function of time in hours (for 0 to 
200 hours) due to the effect of rate constant of Glucose creation, for different values of kg=k=0 .00005 
,0.0005, 0.001, 0.005, 0.05 respectively, (i) Deterministic case (the first upper two row panels), (2) 
Chemical Langevin equation (CLE) (the third and fourth row panels) and (3) Stochastic simulation 
algorithm (SSA) (the fifth and sixth row panels). For CLE we have taken five values of system sizes i.e. 
V=l, 50, 100, 500, 1000; and for SSA we have taken V=l, 25, 50, 75 and 100. 
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Results and discussion 

The dynamical behaviours of p53 and various states of the system induced by glucose concentration 
levels in it are studied using three different computational techniques, namely, deterministic by solving 
equations (2)-(10) where {£} —> 0, CLE by solving equations (2)-(10) and stochastic simulation algorithm 
by simulating the set of reactions listed in Table 2. The parameter values needed for the simulation are 
also given in Table 2. 

Transition of oscillating states: biological rules 

The glucose concentration level (proportional to the rate k which is the rate of creation of glucose given 
in reaction number 9 in Table 2) in the system, due to both synthesized by cellular processes and diffused 
in from extracellular medium, drives the system to different oscillating states (Fig. 2). The p53 dynamics 
maintains stabilized state with minimum p53 level at low values of k (k(0. 00007) and this state is not 
much influenced by intrinsic noise in the system except some random fluctuation about the stabilized 
state (Fig. 2 last four panels of first column). The two dimensional plots of p53 and p53-MDM2 in this 
situation shows fixed point oscillation (oscillation with amplitude zero) (Fig. 4). This stabilized state of 
p53 reflects the normal condition in cellular systems where p53 is maintained minimum level USED. 

The moderate concentration level of glucose in the system (0.0005(fc(0.001) drives the p53 behaviour 
to sustain oscillation state with increasing amplitudes as k increases (Fig. 2 second and third columns). 
In this condition, two dimensional plots of p53 and p53JMDM2 show broaden limit cycle (Fig. 4). 
The random fluctuation in p53 dynamics is minimized in deterministic system and start showing up 
significantly in CLE and SSA where the fluctuations increases as noise strength increases (as V decreases). 
The continuous and periodic changes in the p53 concentration level show the active participation of p53 
in the molecular interaction described in Table 2 which involve creation and degradation of p53 and other 
molecular species in the system. This sustain oscillation state, therefore, may correspond to activated or 
stress state in cellular system induced by glucose. 

The further increase in the glucose concentration level in the system (0.005(fc(0.008) drives p53 dy¬ 
namics from sustain to damped oscillation state and then to stabilized state (Fig. 2 fourth column). In 
this situation, two dimensional plots of p53 and p53-MDM2 show spiral cycle towards fixed point (Fig. 
4). Even though this is a clear damped state in deterministic system and CLE (noise strength is small), in 
SSA (noise strength is comparatively large) the state is still sustain oscilation state with large amplitude. 
This indicates that noise can also be taken as a parameter which can induce stress to the system. 

However, if the glucose concentration level is high enough then the p53 dynamics goes to the stabilized 
state again with high population (Fig. 2 last column). Similary, the comparative behaviour of Glucose, 
ROS, Damaged-DNA and ARF is shown in fig. 3, as a function of time in hours (for 0 to 200 hours) 
due to the effect of rate constant of Glucose creation, for different values of kg=k=0. 00005 ,0.0005, 
0.001, 0.005, 0.05 respectively, (i) Deterministic case (the first upper row panels), (2) Chemical Langevin 
equation (CLE) (the second, third, fourth and fifth row panels) and (3) Stochastic simulation algorithm 
(SSA) (the sixth, seventh, eighth and nineth row panels). For CLE we have taken five values of system 
sizes i.e. V=l, 50, 100, 500, 1000; and for SSA we have taken V=l, 25, 50, 75 and 100. (for different 
colours code see figure 2). The result suggests that as the rate constant of glucose (which corresponds to 
glucose concentration in the system) increaeses the stress within the systems first increases then finally 
acheived an steady state at high value of rate constant. 

In this condition, two dimensional plots of p53 and p53-MDM2 show fixed point oscillation again 
(Fig. 4). In stochastic system, there is still oscillating behaviour at large strength of noise (V = 1) which 
reveals that still noise tries to resist the system to go to stabilized state (Fig. 2 last column fifth panel). 
This second stabilized state may correspond to apoptotic state because increasing glucose level the p53 
dynamics will remain stabilized forever. 

The time of activation T s for p53 temporal dynamics, which can be defined as the time below which the 
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Figure 3. Similar plot as shown in figure 2, of comparative behaviour of Glucose , ROS , 

DamagedoNA and ARF as a function of time in hours (for 0 to 200 hours) due to the effect of rate 
constant of Glucose creation, for different values of kg=k=0. 00005 ,0.0005, 0.001, 0.005, 0.05 
respectively, (i) Deterministic case (the first upper row panels), (2) Chemical Langevin equation (CLE) 
(the second, third, fourth and fifth row panels) and (3) Stochastic simulation algorithm (SSA) (the 
sixth, seventh, eighth and nineth row panels). For CLE we have taken five values of system sizes i.e. 
V=l, 50, 100, 500, 1000; and for SSA we have taken V=l, 25, 50, 75 and 100. (for different colours code 
see figure 2). 


























































Figure 4. Two dimensional plots of p53 for different values of kg in deterministic, CLE (V=100) and 
SSA (V=50). 


dynamics show damped oscillation and above which it shows oscillation death or fixed point oscillation, 
is calculated for different k values (Fig. 6) for deterministic, CLE and SSA. This phase diagram shows 
various regimes of oscillation states and their switching boundaries. The error bars in each curve are due 
to average over 30 ensembles. Further, we calculated amplitudes of the p53 dynamics (A s ) as a function 
of k (Fig. 6 right column) which obeys power law behaviour with k in large k regime i.e. A s (k ) ~ k 1 . 
All the curves of different V are found to be within the error bars showing the similar behaviour of the 
system as a function of V. 

Noise to stabilize the system 

Noise has constructive role in regulating the system by trying to keep the system at normal condition 
(stabilized condition), on removing the noise (as V increases) the system goes to activated or stress state 
(Fig. 5). Further, noise also tries to save the stress state from going to apoptotic state (Fig. 2 fifth 
panels of fourth and fifth columns). It is also evident from the curves in Fig. 6 that as V decreases the 
curves shifts on the right hand side showing (i) helping the system to maintain in normal state by keeping 
away from stress, indicated by larger area occupied towards oscillation death regime as V decreases (left 
hand side of sustain oscillation in Fig. 6 left column) (ii) allowing the system to save from apoptosis by 
slowing down to reach stabilized state (larger area occupied in oscillation death regime in right hand side 
of sustain oscillation in Fig. 6 left column). 

Stability solution 

The stabilized solution of p53 and MDM 2 can be obtained from the stationary conditions of the set of 
equations (2)-(10) by putting -ix = 0 and —> 0,VL Solving for x\ from the nine stationary equations 
we obtain the following equation, 


x 


* 

1 



B 


( 11 ) 
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Figure 5. Two dimensional plots of p53 for different values of V for fixed value of k: (i) CLE for 
V = 1, 50,100, 500,1000 and (ii) SSA for the same values of V. 


where, A = 


/ fcsfcs/cinfei.-ifeii 

(1 1 4 

V ^2^7^12^16 



and B = 


fc4fci2fci6 

2/ciofci3Ali4 ' 


The equation (fill) shows that for large values 


of glucose concentration x 5 (stabilized condition for large glucose concentration), B/x\ —\ 0 which gives 
x* oc y/x|". For small values of glucose concentration, 1 + B/x 5 ~ B/x^, and we found that x\ oc 

maintaining low p53 concentration level. 

Similarly, the stabilized solution for MDM 2 (, x 2 ) is given by, 



x 


* 

2 


c 


B + x % 


( 12 ) 


where, C = y 1 ff) ■ ^ large glucose concentration level, 1 + B/xl —> 1 and equation 

m gives, X 2 oc —4^. Further, for low glucose concentration level 1 + x%/B — > 1 such that xj oc x% 

y x 5 

showing Mdm2 level is at large as compared to higher glucose concentration level. 


Conclusion 

We investigated transition of various oscillatory states in p53-MDM2-Glucose model induced by glucose 
m- Our simulation results in three approaches namely deterministic, CLE and SSA show three distinct 
states, namely, oscillation death, damped and sustained oscillatiory states, and a clear transition among 
these states induced by glucose concentration level. This transition could be the signature of transition 
of p 53 and MDM2 states from normal to stress, stress to apoptotic state induced by glucose. 

Intrinsic noise associated with the system dynamics helps the system to maintain its stabilized state 
and further helps to protect from apoptosis. However, there are other several points to be studied further 
on the roles of noise, finding critical strength of noise below which the role of noise is constructive, 
role of external noise in maintaining orderness in non-equilibrium and non-linear physical, chemical and 
biological systems. 
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Figure 6. The phase diagram for different V in ( T s — k) plane and {A a — k) are presented for (i) 
deterministic (V oo), (ii) CLE (V=l, 25, 50, 75, 100) and (iii) SSA (V=l, 25, 50, 75, 100). 
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